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ABSTRACT 

We show that the MBH-A^buigc scahng relations observed from the local to the high-z Universe 
can be largely or even entirely explained by a non-causal origin, i.e. they do not imply the need 
for any physically coupled growth of black hole and bulge mass, for example through feedback by 
active galactic nuclei (AGN). Provided some physics for the absolute normalisation, the creation 
of the scaling relations can be fully explained by the hierarchical assembly of black hole and stellar 
mass through galaxy merging, from an initially uncorrelated distribution of BH and stellar masses in 
the early Universe. We show this with a suite of dark matter halo merger trees for which we make 
assumptions about (uncorrelated) black hole and stellar mass values at early cosmic times. We then 
follow the halos in the presence of global star formation and black hole accretion recipes that (i) 
work without any coupling of the two properties per individual galaxy and (ii) correctly reproduce 
the observed star formation and black hole accretion rate density in the Universe. With disk-to- 
bulge conversion in mergers included, our simulations even create the observed slope of '^1.1 for the 
-^'^BH^-^'^buigc-relation at z = 0. This also implies that AGN feedback is not a required (though still 
a possible) ingredient in galaxy evolution. In light of this, other mechanisms that can be invoked to 
truncate star formation in massive galaxies are equally justified. 

Subject headings: galaxies: fundamental parameters — galaxies: nuclei — galaxies: bulges — galaxies: 
evolution 



1. INTRODUCTION 

About a decade ago tight correlations between 
galaxy properties and those of central supermas- 
sive black holes (BHs) were empirically established: 
BH masses scale with the luminosity, the mass 
and the velocity dispersion of their host galaxies' 
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quickly interpreted to yield two important implications: 
(i) If these correlations were to exist for a random 
set of galaxies, then every galaxy should conta in a 
supermassive BH fe.g. lKormendv fc Richstonelll995| ) (ii) 
If every galaxy contains a BH and given the observed 
scaling relations, the evolution of galactic bulges and 
central BHs should be coupled by a physical mechanism 
(the 'co-evolution' picture). 

The former statement appears to be true at least above 
some mass threshold, and introduced a new ingredient, 
the central BH, in galaxy evolution. The latter state- 
ment is based on the vast energy available from accreting 
black holes, providing the easiest conceivable mechanism 
to physically couple BH and bulge properties despite the 
difference in linear scales. Coupling only a few percen t of 
this energy in an "AGN feedback" (ISilk fc Reesl(T998h to 
the gas within the galaxy would have vast implications 
on the temperature and structure o f the surrounding in- 
terstellar medium. Ad hoc models (Graiiato et al.l 120041 : 
iDi Matteo erall[2005l : ICroton et aimo^ were very suc- 
cessful in generating feedback loops that involve the en- 
ergy from AGNs for quenching star formation (SF) and 
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fueling of the AGN themselves. Different incarnations of 
AGN feedback are in principle able to not only couple 
SF and BH accretion, but to simultaneously fix a num- 
ber of existing problems in galaxy evolution, namely an 
overproduction of massive galaxies in semi-analytic mod- 
els as well as the inability to truncate SF fast enough 
to reproduc e the observed col or-magnitude bimodality 
of galaxies (jBaldrv et al.ll200^ . This motivated the in- 
clusion of AGN feedback in a yet-to-be-determined form 
and physical description as the driving force behind the 
BH-bulge scaling relations. 

Although the actual effectiveness and impact of at 
least "quasar mode" feedback models is still unclear, 
the interpretation of the scaling relations as a physi- 
cally coupled evolution is largely assumed to be cor- 
rect and continues to provide the basis for many stud- 
ies. As an examp l e, the evolution of th e scaling relations 
dTreu et al.' 12004; 'Peng et al.' '2006a'.'b'; 'Treu et al.l 120071: 
i Schramm et al. 2008; Jahnke et al. 2009; Merlon i et al l 
I2OIOI: iBennerte't al.ll2010l: iDecarU et al.ll20ia) is investi- 
gated in order to constrain the physical drivers behind 
co-evolution and the growth mechanisms of BHs. 

1.1. An alternative origin of the scaling relations 

IPend (|2007[ ) demonstrated a thought experiment for 
the potential origin of a MBH-.Wbuigc-i'elation without a 
physical coupling, but as the result of a statistical con- 
vergence process. In short, he showed that in principle, 
arbitrary distributions of MBH/Af*-ratios in the early 
universe converge towards a linear relation through the 
process of galaxy merging. In this central-limit-theorem 
view, a large number of mergers will average out the ex- 
treme values of Mbh/-^* towards the ensemble average. 
What was deliberately left open in this experiment was 
whether there are enough major galaxy mergers in the 
history of an actual galaxy in order to drive this process 
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far enough. 

In this present study we pick up this thought by fol- 
lowing realistic ensembles of dark matter halos through 
cosmic time. Our immediate aim is to test whether the 
simple assembly of galaxies and their BHs according to 
a ACDM merger tree is able to produce BH-galaxy scal- 
ing relations from initial conditions at early times, where 
Mbh and Mbuige (or M*) were completely uncorrected 
per individual galaxy. We circumvent the inherent prob- 
lems and degrees of freedom of a full semi-analytic model 
by not trying to simultaneously solve the problem of SF 
truncation or correct BH or stellar mass function, but 
restrict the question solely to the genesis of the scaling 
relations. As an input for SF and BH accretion we use 
observed relations, but prevent any recipes that im- or 
explicitly couple BH and stellar mass growth per individ- 
ual galaxy. 

2. NUMERICAL SIMULATIONS AND MERGER TREES 

In this paper we u se the Lagrangian code PINOCCHIO 
([Monaco et al.ll2002[ ) to construct high resolution ACDM 
merger trees (f or a comparison between Nbody codes and 
PINOCCHIO seelLj et al. 2007). The simulation has a box 
size of 100 Mpc, and 1500'^ particles, this ensures a very 
high mass resolution, rup = 1.01 x 10^ Mq. The construc- 
tion of merger trees is straightforward with PINOCCHIO: 
the code outputs a halo mass every time a merger occurs, 
i.e., when a halo with more than 10 particles merges with 
another halo. From an initial 6.5 x 10® halos, we receive 
a resulting sample of 10932 halos with M > 10" Mq at 
z = 0. 

When two halos merge, the less massive one can ei- 
ther survive and continue to orbit within the potential 
well of the larger halo until z = 0, or merge with the 
central object. The averaging process described above 
will only apply to this second category of halos (the ones 
that actually merge), we then adopt the dynamical f ric- 
tion formula presented in lBovlan-Kolchin et al.l ()2008[ ) to 
compute the fate of a halo. The orbital parameters of 
the halo are extracted from suitable distributions that 
re produce the results o f Nbody simulations as described 
in [Monaco et al.l (|2007f ). If the dynamical friction time 
is less that the Hubble time at that redshift, we consider 
the the halo to merger at a time t = tdym if it is longer 
the satellite halo is removed from our catalog. Each halo 
in our sample at z = has formed by at least 200 mergers 
and the most massive ones have had more than 5 x 10** 
encounters. 

3. CREATING SCALING RELATIONS: AVERAGING AND 
MASS FUNCTION BUILD UP BY HIERARCHICAL 
MERGING 

The main message of this work is to demonstrate which 
effect merging over cosmic time has on an ensemble of 
halos with initially uncorrelated A/bh and values. 
These values change their distribution and converge to- 
wards a linear relation by z = - in the absence of SF, 
BH accretion, disk-to-bulge conversion and hence any 
physical connection between the two masses. 

For this task we follow dark matter halos through their 
assembly chain. We assign a stellar and a BH mass 
to each dark matter halo once its mass becomes larger 
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Fig. 1. — Changes of Mbh vs. M* from an initially uncorre- 
lated (within 4 dex in each parameter) distribution at high z (blue 
points) to z = purely by mass assembly along the merger trees, 
i.e. without SF, BH accretion and disk-to-bulge mass conversion. 
A very tight correlation of slope 1.0 is created by the merging, 
with smaller scatter for higher masses which experienced more 
merg er. The black line is the observed local AfBH-A^*-relation 
from lHaring fc Rixl 12001 ') with slope=1.12. 

than 10* M(^, the corresponding redshift in the follow- 
ing is called Zf. We set our initial guesses for and 
Mbh as a fixed fraction of the dark matter mass plus a 
(large) random scatter. We used M*/Mdm = 10^"^ and 
Mbh/A^dm — 10^*" for the initial ratio; the scatter is 
taken from a logarithmically flat distribution of 3 dex 
for the two quantities (blue squares in Figs.[T] and[2|). 

We have no knowledge of any realistic seed mass scat- 
ter, but take 4 orders of magnitude variations in the 
A/BH/Af*-ratio as a proxy for "uncorrelated". Em- 
pirical constraints on the possible parameter space for 
seed black hole mass d o not seem to support s eeds more 
massive than ~10^ M^ HVolonteri fc Nataraianl2009) . to- 
wards lower masses a few solar mass BHs are clearly be- 
ing produced by stars. Whether this matches the true 
distribution of seed masses is not important, but by sim- 
ply taking a large range that is currently not ruled out 
represents a rather conservative starting point for our 
demonstration. Halos are then propagated along the 
merger tree to z = (the red points in Figs. [T] and 
[2]). When two halos merge according to our dynami- 
cal friction formula, we set the resulting stellar and BH 
masses equal to the sum of the in dividual masses before 
the merger (|Volonteri et al.ll2OO30 . The final mass m 
A/bh and as well as the corresponding normalization 
is determined simply by the sum of the individual halos 
contributing to a final halo. 

Fig. [1] shows that the hierarchical formation of galaxies 

^ We picked this m ass since at lower m asses halos likely did/do 
not form stars at all lIMaccio et al.|[2010l) . The most massive pro- 
genitor of 2; = galaxies form according to this definition in the 
range z = 15 — 17, while low mass satellites can form as late as 
z ~ 3. 
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provides a strong inherent driver from the uncorrelated 
initial distribution to a linear relatior0. This effect is 
independent of the chosen initial conditions as Figure [2] 
demonstrates, where completely different initial condi- 
tions result in a relation with the same slope and very 
similar scatter. 

This experiment shows that the dominating structural 
parts of the observed MBH-Af*-scaling relation - i.e. (1) 
the existence of such a correlation, (2) that it extends 
over several orders of magnitude in mass, (3) the fact 
that the slope is near unity, and (4) an increasing scat- 
ter to lower masses - can be explained by this physics-, 
feedback- and coupling-free process. A slope^l scaling 
relation does not need any physical interaction of galaxy 
and black hole. In the next sections we will show that 
this holds also when adding "2nd order" effects like ac- 
tual star formation and black hole accretion, as well as 
disk-bulge conversion. 

4. ADDING STAR FORMATION, BLACK HOLE 
ACCRETION, AND DISK-TO-BULGE CONVERSION 

We so far demonstrated that merging alone is the basic 
mechanism to create a Men-Af* scaling relation. How- 
ever we need to add a number of ingredients to our model: 
Placing all mass already at high redshift is not conserva- 
tive, since all of stellar and BH mass is subject to the full 
merger averaging process. In the actual universe we know 
that the majority of BH and stellar mass in the universe 
was c reated after z ~ 6 (jSoltanl ll982': 'Hop kins fc Beacoml 
|2006[) . and hence experiences less merger generations. 
Moreover pure merging produces a monotonic relation 
between Mdm and M^, as shown in the upper panel of 
Figure [21 which is at odds from empirical results. Since 
SF and BH accretion density depend on redshift, we will 
add the approximate right amount of SF and BH growth 
at the right cosmic times, and thus at the right place in 
time with respect to the merger cascade. The goal of 
this excercise is not to create a full semi-analytic model 
of galaxy formation, but to test which effect realistic 
assumptions about mass growth have on the resulting 
A^BH-Af* scaling relation. 

To construct SF and BH accretion recipes we will use 
three observed relations as input: (i) The halo occu- 
pation distribution, i.e. the rel ation between dark mat- 
ter halo mass A/dm and (jMoster et al.l i2010'). (ii) 
the Lilly-Madau relation for the evolution of SF rate 
(jHopkins fc Beacoml I2006D . an d (iii) the evolutio n of 
bolometric luminosity of AGN ([Hopkins et al.ll20Q7[ ). 

We want to explicitly note that by using these rela- 
tions, we do not force any coupling of Mbh and or 
Afbuigc pel" individual galaxy, but all recipes only relate 
to ensemble averages. Any potential implicit couplings 
only act on the ensemble and not on an individual galaxy, 
hence they can not induce a correlation in the (originally 
uncorrelated) data points. 

4.1. Star Formation 

^ The convergence is in fact too strong (see next sections), as 
the scaling relation it produces by 2 = is much tighter than the 
observed 0.3 dex scatter. In principle the scatter has a depen- 
dency on the number N of merger generations, but the relation gets 
complicated by the different masses of the merging components and 
different merging times across the tree. 



Up to now, at redshift z = the stellar mass of our 
galaxies is simply given by the sum of all stellar masses 
of its j-progenitors M°(i) = ^M^{zf){j), masses that 
as explained were originally drawn from a random dis- 
tribution. As shown in the upper panel of Fig. [31 for 
a given halo mass the stellar mass obtained in this way 
is too low when compared to the empirical expectations 
from Halo Occupa tion Distribution (HOD) models (e.g. 
iMoster et al.ll2010() . This is because the total galaxy is 
more than the sum of its seeds and we have neglegted 
star formation so far. On the other hand. Fig. [3] tells 
us exactly how much stellar mass each halo is missing. 
Hence we take this constraint, the HOD results, to fix 
the stellar mass produced throught star formation: 

M*sp{i) = M*{Mdmii)) - M*{i) (1) 

where M* {Mdm(i)) is the expected stellar mass for the 
z-th halo with dark matter m ass Mdmii) as predi cted by 
the HOD model presented in lMoster et all (|2010D . 

Now we need to distribute this stellar mass from SF 
along time, i.e. among all progenitors of galaxy i along 
the merger tree. We do that according to the following 
formula that gives the stellar mass produced through SF 
for the j-th halo in the merger tree of final i halo: 

M/^ = ^ X M^Jj) X LT{]) X /(z/(j), zm(j)). (2) 

The constant A is fixed by the requirement that 
M^^{i) = and is the same for all progenitors. 

The time LT is the lifetime of a halo, defined as time 
between the formation redshift {zf. when A/dm > 10® 
Mq) and the moment Zm, when it merges with a more 
massive halo, which is not necessarily the main branch 
of the merger tree. With this definition we assume that 
galaxies are able to actively form stars only when they 
are the central object within their host halo. 

The function / is used to give different weights to the 
life time LT at high ad low redshift, in this way, for a 
given life time, a galaxy will produce more (less) stars at 
high (low) redshift, according to the Lilly-Madau plolQ. 
We define f{zi,Z2) as the integral between zi and Z2 of 
the assumed star formation rate (SFR): 

fizi,Z2) = f 'sFR(z)dz. (3) 

In our reference model we assume for the red- 
shift evolu tion of the SFR the func tional form sug- 
gested by [Hopkins fc Beacoml ([2006[ ). namely results 
listed in Table [I] for a modified Salpeter IMF (see 
lHopkins~fc Bcacom 2006, for more details). Finally the 
factor AT* (z) takes care of the observed mass-dependence 
of specific star f ormation rates and we fixed the expo- 
nent q = 0.8 (e.g. lDaddi et al.[[20Q7[: [Bouche et alB OOQ). 
In the Appendix [^ we will present results for different 
choices for q and SFR{z). 

Let us summerize one more time our parametrization 
for star formation: When the j halo appears at z/(j) it 
gets an initial stellar mass {M^,{zf){j)) from a random 
distribution as described in Section [3| Then it will "pro- 
duce" its own stars (Mj^) until it will be accreted onto, 

* We note that we do not include a different shape of this star 
formation history as a function of mass. As the Appendix shows, 
this will not have any effect on the results 
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Fig. 2. — The scaling relations are produced independently of the initial conditions: Shown are four vastly different initial conditions 
— for a subset of our halos for better visibility all leading to a slope=l relation at 2 = with similar scatter. The 2 = normalisation 
comes out differently since the geometric mean of the initial masses is different (flat distribution in logarith for both quantities). Symbols 
and the line have the meaning as in Figure [T] 



and become part of, a more massive halo. During its life- 
time it will also accrete stellar mass from merging (lower) 
mass haloes. If halo j will merge with the central galaxy 
it will add a fraction of its stellar mass to the bulge of 
the central galaxy as described below in Section If 
halo j will merge with a more massive halo (k) before 
merging with the central halo it will give all its stellar 
mass to halo k and cease to exist as a halo of its own. 

4.2. Disk to bulge conversion 

We assume that all stellar mass produced through star 
formation will occur in the disk component of each halo 
and then apply a recipe to convert part of this disk mass 
to bulge mass as a consequence of mergers. 

The amount of disk-to-bulge mass conversion depends 
on a multitude of parameters as mass ratio, gas frac- 
tion, and orbital parameters of the merger, which are 
impossible to implement in our context. Instead we fol- 
lo w a simpler recipe inspired by the numerical results 
of iHopkins et al.. (,2009. ) . solely depending on the stellar 



mass ratio of the two merging partners: (i) Bulge mass of 
main halo and satellite will just be co-added, (ii) the disk 
mass of the satellite fully goes into the resulting bulge, 
and (iii) a fraction of the main halo disk, directly pro- 
portional to the mass ratio, also gets converted into bulge 
mass. With this recipe we are able to approximately re- 
produce the ratio of bulge to total mass observed in the 
local universe. 

4.3. Black hole accretion 

Since the mechanisms of BH accretion continue to be 
unclear, we assume a simple recipe; the BH will double its 
mass in a stochastic way on a characteristic time scale 
T. This will happen for all j progenitors of halo i. 
Similar to what we have done for star formation, we link 
the number of doublings of a given halo to its life time 
and we weight this time with a function g similar to the 
function / in Section |4T] 

In practice the number of mass doublings of the j-th 
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Fig. 3. — The relation between M* and MpM at z = 0. The 
black line shows the predictions from lMoster et aLl 12010) . the red 
point the results from our simulations. Upper: Only merging but 
no SF taken into account. Lower: With our SF recipe implemented. 

black hole in the i merger tree is given by the expression: 

iVdoubj - LT{j) X g{zf{j), Zrn{j))/T, (4) 

Where LT has the same meaning as in Eq. [2l For 
black hole accretion we choose as the weighting function 
g(zi : Z2) the integral between zi and Z2 of the bolomet- 
ric luminosity of AGN {AGNl{z)): 



9{ZI,Z2) 



AGNL(z)dz. 



(5) 



In our reference model the redshift evolution of 
the AGN bolometric luminosity is modeled with a 
double power law a imed to reproduce the results of 
iHopkins et all (|2007l . Fi gure o : 

( AGNl{z) \ _ \ 2.02 ■ log(z) + 7.83 for z < 1.7 
UoMpc-^ y I -2-09 • log(z) + 8.78 for z > 1.7 

Similarly to the / function, g can be used to allow for 
higher accretion rates at high redshift compared to low 
redshift for a fixed halo life time. 

The characteristic time t is chosen in order to match 
the normalization of the observed AfBH~-Mbuigo-scaling 
relation at z = for log(MBH) = 7 and we obtained 
r = 1.9 X 10^ Gyrs. As a stochastic element to BH growth 
we cast for each of the N^ouh.j events a random number 
from a flat distributing in the range [0:1] and effectively 
double the BH mass only if this random number is > 0.5. 
For the main branch the number of doublings is of the 
order of 6-8, while it is in the range 0-3 in the other 
branches of the tree. In the Appendix El we will present 
results for different choices of the parameters involved in 
equation ID 

4.4. Resulting scaling relations 



The Mbuigo~AfBH-distribution at z = with the above 
receipes added is shown in Figure [H with the observed 
values overplotted. Compared to the pure merger assem- 
bly in Figures [1] and [2] we note a substantially increased 
scatter - closer to the observed - and a steeper than lin- 
ear relation. Still, despite both the shift of SF and BH 
growth to later times as well as the random parts of SF 
and BH accretion, a clear correlation is produced. The 
effects of SF, BH accretion and disk-to-bulge conversion 
do not destroy the correlations but only induce a "2nd 
order" modification of them. Over time a mass function 
is built up and the scatter decreases substantially from 
the initial 4 dex, particularly for the log(Afbuigc) > 10 
regime. 

The simulated relation is very similar to the observed 
points, it reproduces the slope >1 almost perfectly, even 
the fact that at the high mass end the observed points lie 
above the mean slope - and this without adjustable pa- 
rameters beyond normalization. Also, the higher scat- 
ter at low masses can be seen in both simulations and 
observations, the "cloud" above the z = relation in 
our model data represents disk galaxies with small bulge 
masses as obs erved by Greene ct al. (2008). As can be 
seen in Figure IATH in the Appendix, these objects would 
still lie on the local mass scaling relation, but with their 
total and not their bulge mass. 

We want to stress that the results in Figure |4] do not 
depend on our parametrization of SF rate or BH accre- 
tion. For example changing the functional form of / or 
g in Eqs. [2] and |4] only marginally affects the scatter of 
the simulated MeH-Afbuige-relation and leaves the slope 
unchanged. The same is true for the other parameters 
described in the previous sections, see the Appendix for 
details. There is a sole exception to this, the assumed 
initial seeding masses: The larger the initial A/bh and 
Af*, the less mass has to be created by SF and BH ac- 
cretion. In this way less mass is entering the halos at 
later times and more mass is subject to the full cascade 
of mergers, which will lead to a smaller scatter in the 
scaling relations at z = 0. 

5. DISCUSSION 

We showed above that all basic properties of the BH- 
bulge mass scaling relation in the local Universe - a rela- 
tion between properties of individual galaxies - are pro- 
duced naturally by the merger-driven assembly of bulge 
and BH mass, and without any coupling of SF and BH 
mass growth per individual galaxy. 

The convergence power of galaxy merging is very 
strong for a realistic halo merger history, even with a 
correct placement of SF and BH accretio n along cos - 
mic time. This means that the mechanism iPend ()2007f ) 
sketched in his thought experiment works also in a real- 
istic Universe - there is enough merging occurring in the 
universe and hence (most of) the scaling relations can 
be entirely explained without any physical mechanisms 
that directly couples Afbuigo and A/bh growth for a given 
object. 

5.1. Implicit coupling of black hole and galaxy? 

The scaling relations are produced naturally in our toy 
model - but does this preclude any implicit coupling of 
Mbh and AI^ by design? The most obvious features to be 
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Fig. 4. — Left: Initial uncorrelated high redshift seeds for M, and A/bh (blue filled squares) and resulting 2 = ^'I'SM^^lhvA^a scaling 
relation for a subset of 400 randomly selected merger trees (red points), compared to the observed local relation in black, including the 
com pilation from (Feoli Mancini 2003 . circles), low- mass spheroids (open squares) and upper limits for spiral bulges (triangles) from 
iGreene et al.i (|2008|). The solid line is the linear fit by IHaring fc Rixl (|200D with a slope of 1.12. Right: The full set of resulting 10932 
galaxies at z = with the low— z data overplotted. 



discussed in this respect are (a) the shape of the halo oc- 
cupation distribution and (b) the question, which mech- 
anism determines the relative value of BH accretion to 
SF, hence the absolute normalisation of the MBH-.Mbuigo- 
relation. 

(a) HOD shape: The HOD was empirically inferred 
and shows that the ratio of stellar to dark matter mass 
is not constant but is a function of mass itself. Towards 
the massive end stellar mass does not increase in parallel 
to DM, star formation appears suppressed. Its impact on 
the MBH~-^^buigc-relation is the slight curvature in Fig- 
ure 21 with the noticable upturn at the massive end - 
consistent with the observations. The HOD shape rep- 
resents a long known feature of galaxy formation and 
was "fixed" in models by introduction of a quenching 
mechanism, suppressing SF above some mass, often by 
inclusion of AGN feedback receipes. 

In our toy model however, we do not make specific as- 
sumptions about which mechanism produces the HOD 
we use as a constraint. We do not tie BH accretion or a 
merger event to the suppression of SF - it can be sup- 
pressed by any mechanism, e.g. a modified SN feedback 
receipe or gravitational heating by infalling clumps of 
matter. The latter mechanisms are completely indepen- 
dent of BH accretion and, while modifying the HOD, 
by nature can not have an impact on the MBH~-Mbuigc- 
relation. Even if AGN feedback was the source for shap- 
ing the HOD, this would only be the cause of the second 
order shape deviation from a linear slope, not the exis- 
tence of the MBH--Mbuige-relation itself. 

(b) Absolute normalisation: Our initial model just 
propagates the (rather ad hoc) seed masses in stars and 
BH to z = 0, i.e. the normalisation of the MBH--Wbuigo- 
relation at z = is directly set by the ensemble mean 
ratio of seed masses. Since most of stellar and BH mass 



in reality is produced by SF and BH accretion later on, 
the high-z ratio is in fact unconnected to the z = nor- 
malisation. 

In our simulation we set the normalisation by request- 
ing a match of our simulation results with the empir- 
ical MBH-Afbuigo-relation at M*=1O^°M0. Arguments 
have been brought forward that the actual normalisation 
must be the result of a regulatory feeback loop involving 
BH accretion and SF. This is an attractive scenario, as 
it would explain both the creation of the scaling rela- 
tions as well as their normalisation. Models of this kind 
have been implemented in several semianalytic models of 
galaxy formation, in all cases with free parameters that 
actually control the absolute normalisation of the result- 
ing scaling relations, set to ad hoc values to again match 
this and other observations. Since we show in this paper 
that (most aspects of) the AfBH^-^buigc-scaling relations 
are created automatically by hierarchical assembly, these 
feedback models actually seem to achieve too much - the 
creation of a certain Mbh/-^* ratio for each galaxy and 
at all times, which as a conspiracy would come on top of 
the formation path of the scaling relations demonstrated 
here. 

This said, we want to sketch the outline of an alter- 
native scenario, that could well be responsible for the 
absolute normalisation but has not yet been explored: 
The main ingredients for both SF and BH accretion are 
(i) gas, and (ii) a trigger to form stars or to bring down 
gas to the BH. For both stars and BHs the amount of 
growth is basically a product of the two ingredients. 
At early times gas was ample and the number of both 
galaxy mergers and gas disk instabilities was high un- 
til the peak of activity around z = 2. The triggering 
mechanisms subsequently decreased with the decreasing 
number of mergers and reduced gas reservoirs, either in 
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number or duration or both. If both SF and BH accre- 
tion were to be ruled by a set of random triggering mech- 
anisms and the specific gas fraction in a galaxy, then very 
high-z galaxies might exhibit a strong variance in their 
-^^BH/Af*-distribution as well as in their (instantaneous 
and also time averaged) BH accretion over SF ratio, with 
a dependency on the actual total mass, morphology, gas 
mass, trigger type, environment, etc. However, as a cos- 
mic mean there will be a global Mbh/M^ value, just as 
an average over the efhciency of the spectrum of ran- 
dom triggering mechanisms and realised conditions to 
produce new stars or BH mass - a number which could 
also change with time. 

The hierarchical assembly of galaxies and its averag- 
ing mechanism now relieves us from having to search for 
regulatory mechanism per galaxy, since with each galaxy 
merger the spectrum of Mbh/M^ values will increasingly 
average out. The fact that the slowly starting deple- 
tion of gas reservoirs in galaxies at z ~ 2 is accom- 
panied with, and not preceded by, a slowly decreasing 
merger rate has the consequence that at higher redshifts 
there were enough mergers to average out the extreme 
AIqh/M^ systems - at lower z, with a decaying gas reser- 
voir and merger rate and the transitio n to a "secular uni- 
verse" (see e.g. iCisternas et aLll2011[ ). the prerequisites 
for producing new extreme values become less and less 
frequently fulfilled. The ensemble converges towards the 
observed MBH-A<rbuigo-relation at z = 0. 

Recently, first nested multi-scale simulations of gas in- 
flow into the very centers of galaxies have been successful 
(|Hopkins fc Quataertll2010t ) and give a first impression of 
how in principle random instabilities, strongly depending 
on the actual conditions in the galaxy, can create gas in- 
flow into the very center and the fuelling of either BH, 
SF or both. It is mechanisms of this kind that will cre- 
ate a certain mean value of Mbh/M^, averaged over all 
galaxies, which can be vastly different for an individual 
galaxy. How different is so far unclear, and whether this 
mechanism for BH fuelling precludes "runaway" growth 
of BHs to 100- or 1000-fold in a single instance, though 
unlikely, needs to be seen. All of this can in principle 
be realized without any AGN feedback - though it does 
not rule it out -, and has the freedom to have a mass- 
or environmental-dependent efficiency component, as we 
also observe a non-unity slope of the z = scaling rela- 
tion. 

Coming back to the initial question, we conclude that 
there is no necessity that our toy model makes any im- 
plicit assumptions about a physical Mb h~-^^*- coupling. 

5.2. Further consequences of this mechanism 

Our mechanism is also able to explain other observa- 
tional results that are often used as evidence to support 
the picture of agn feedback. 

(a) One of them is the potentially different Mbh- 
Mbuige-re lations in gal ax ies w ith classical and pseudo- 
bulge s (|Greene et all 120081: iGadotti fc KauffmannI 
iOOgl). Pseudo-bulged are thought to be formed 
through secular processes, ra ther than major merging 
(jKormendy fc Kennicuttll2004f ). This has the immediate 
implication that the bulge has taken a different long- or 
mid-term assembly route compared to the BH, hence it 
is actually not expected that galaxies with pseudo-bulges 
obey the same MBH-Mbuige-relation than those with 



classi cal ones. 

(b) iHopkins et al.l ()2009D found that the stellar mass 
inside a very small radius near the black hole M*(< R), 
comparable to the BH's sphere of influence, shows a much 
larger scatter than the scatter in AIbk/M^,. They argue 
that this is an indication that gas was transported to near 
the BH, to form said stars, but that this had apparently 
no impact on the scatter in Mbh , hence a self regulating 
mechanism should have been at work. 

With our model also this can be simply explained: The 
BH and bulge of a galaxy take part in the same long-term 
assembly and averaging cascade, hence their values cor- 
relate and scatter in their ratio is small. The central 
density of stars on the other hand does not. It is gov- 
erned by more short term gas inflow, star formation and 
redistribution mechanisms, hence it is expected that it 
does not correlate well with the overall mass of the bulge 
or the BH. 

(c) Following this line of argument, we in princple 
expect a correlation with Mbh for any (mass) param- 
eter that is subject to the same ACDM assembly chain. 
This includes the halo mass, to some extent the total 
ma ss of the galaxy (for both see Appendix [B] and Fig- 
ure [A3|, but also e.g. the total mass of globular clusters 
in a galaxy, which recently ha s been found empirically 
bv lBurkert fc Tremaind (|2010D . 

(d) On the other hand, our model is too basic to 
be able to re produce higher order effects. A number 
of stud ies (AUcr fc Richstonc 2007: Hopk ins et al.l [20071 : 
iFeoh fc Mancini 2009; F eoli et al . 2010) have suggested 
that the most fundamental relation with Mbh is neither 
A^buigo nor (Jbuigo but that rather galaxy binding energy 
or potential well depth. 

What these studies actually find are residual correla- 
tions in the scaling relations that are in some way related 
to the the compactness of a bulge/spheroid at a given 
Mbh- This can be effective radius or binding energy or 
any other quantity that is ex- or implicitly a measure of 
compactness. Since bulge mass depends just linearly on 
the progenitor mass ratios, the exact size, compactness, 
binding energy of a spheroid produced in a merger will 
depend non-linearly on other parameters like gas frac- 
tion, merger orbit, etc. These parameters could easily 
be responsible for the residual correlation found in the 
MBH~-^buigo-relation as they are measures of the specific 
short- or mid-term merger history of each galaxy, while 
Mbuigo is a long term integral. Our deliberately simple 
toy model itself can not make any statements on this 
issue as it does not trace e.g. galaxy scale radii. 

Interpreting the above points with respect to the rele- 
vance of AGN feedback, we find no strong argument for 
AGN feedback as a necessary mechanism at work. How- 
ever this does not mean that AGN feedback does not 
exist. It only means that AGN feedback is still a possible 
mechanism involved in parts of galaxy evolution, which 
might be important e.g. for subclasses of the galaxy pop- 
ulatior0. This implies that other proposed means of en- 
ergy (or momentum) injection that h ave been proposed 
to quench SF in massive galaxies fe.g. lDekel fc BirnboimI 

^ Undoubtedly, at least "radio mode" feedback l|Croton et al.l 
200611 is actually observed in s ome massive clusters (e.g. 
Fabian e t al. '200ff; 'Be st et al.|[200a 'l and possibly on the group- 
level l|Gio dini et al, 20l0). 
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[200l:IKhochfar fc Ostrikeil2008t EoFaro et alJl2009D are 
viable options. In other words, all other mechanisms 
that can be invoked to truncate star formation in mas- 
sive galaxies appear to be equally justified. 



5.3. The comparison with previous work 

Our result is different to all previous studies, and with 
more than a decade of semianalytic models including BHs 
having passed, one of the obvious questions is: Why 
have others not found this before? The closest studies 
to our s is a short proceedings contribution bv iGaskelll 
(|201C1| ). whic h is basically paraphrasin g iPend ()2007[) . and 
the work bv iHirschmann et al.l (|201Clf ). Their study is a 
very systematic assessment of how the scaling relation 
scatter evolves under the influence of galaxy merging. 
However, their initial setup in all cases was an existing 
correlation of Mbh and , and due to their specific goal 
they deliberately ignored the influence from SF and BH 
accretion. 

Other recent studies investigating thi s subject ei- 
ther e xplicitly include AGN fe e dback Ce.g. iCroton et al.l 
200l iRobertson et all [200l iBooth fc Schavd 120091: 
Johansson et al.l 120091 : IShankar et alJ l2009t ) or couple 



5.4. Impact on evolution studies 

In the last decade a substantial number of at 
tempts were made to measure local and higher red 
shift scaling relations of MBH~-^^ buiirei AfsH-Cbuige . 
or MsH-I ' buire, (iTreu et al. 2004; Peng et al.l l2006allb| : 
Woo et al.ll2006t iTreu et al..i2007[ iSchramm et al.l 1200 



the merging scenario wi th a self regulation prescrip- 
tion for BH growth (e g. iKauffmann fc Haehneltl 120001: 
iVolonteri fc NataraianI 120091. In retrospect it is quite 
obvious why they did not find our results earlier, as 
most of these models were developed for galaxy evolu- 
tion in general. Once BHs entered the equations, they 
added terms of (regulated) BH growth and, when eval- 
uating the generated MBH^.^^buigc-relations, noted that 
their prescription managed to produce a decent match 
to the empirical relation. The interpretation that this 
meant the AGN feedback- or coupling-prescription were 
correct, now requires re-evaluation from our new point of 
view - it was the simple result of hierarchical assembly 
at work. 



Jahnke et all 120091 iMerloni et al.l 120101 : iBennert et al l 
20101 : iDecarli et al.l '2010*) and to interprete them with 
respect to the mechanisms that couple BH growth and 
their impact on galaxy formation. Which implications 
does the non-causal origin of the scaling relations have 
for these results? If the mean cosmic M* and Mbh actu- 
ally evolved similarly, our results explain at least a part 
of the bulge scaling relation evolution for galaxies with 
substantial disk components: it is the simple conversion 
of disk to bulge mass in galaxy mergers. What still re- 
mains interesting and needs to be substantiated is how 
at high redshifts the relation between Mbh and total M* 
(or e ven Mhnire for bu lge dominated galaxies) evolves 
(e.g. IWalter et al.ll200^ . This would continue to pre- 
dict a substantial early BH growth - with corresponding 
implications for BH feeding models. 

One other aspect that could serve as a diagnostic is the 
evolution of the scaling relation scatter. When coupled 
with predictions of BH and stellar mass assembly from a 
proper model, the scatter can be used to study e.g. the 
distribution of seed Mbh at early times. We will follow 
up on this issue in a future publication. 
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APPENDIX 

TESTING THE MODEL: ARE THE PARAMETER CHOICES SPECIAL? 

Few (free) parameters enter in our parametrization of star formation, black hole accretion and bulge to disk conver- 
sion. In this section we want to explore different choices with respect to our reference model and test their impact on 
the final resu lts. All different models are listed in Table [T] 

Figure \KA\ shows the effect of varying our 'weighting' functions. In model B we set SFR{z) = 1.0, e.g. we assume 
a constant star formation rate as a function of redshift; in model C we set also AGNl{z) — 1.0. The resulting A^bh- 
Mbuigc-relations are indistinguishable from the original (A) model. The bottom right panel of figure lA.ll (model D) 
shows an even tighter correlation between AfBH~-^^buigcwith respect to our reference model. This is because in model 
D we remove the stochasticity in the BH mass doubling, forcing all BHs to double their mass every r Gyrs. This 
increases the number of doublings in the merger tree branches with short life time, increasing the fraction of BH mass 
that is accreted through mergers and hence is subject to the central limit theorem. 

In Figure|X2]we check the effect of our parametrization of dynamical friction (model E, where the dynamical friction 
time is multiplied by five) and of our disk-to-bulge conversion, assuming that a fraction of the disk mass proportional 
to the square- root (model F) or proportional to the square of the merger ratio is promoted into the disk (model G). 
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TABLE 1 

Model parameters 



Modcl^ 


SFR{z)b 




AGN£{z)d 


Stochasticity° 


B/D conversion' 


Dynamical Friction^ 


A 


HB06 


0.8 


H07 


Y 


Mgi/Mg2 


BOB 


■R 
JJ 


const. 


0.8 


H07 


Y 


M-G1/MG2 


JjUo 


C 


const. 


0.8 


const. 


Y 


Mgi/Mg2 


BOB 


D 


HB06 


0.8 


H07 


N 


Mgi/Mg2 


BOB 


E 


HB06 


0.8 


H07 


Y 


Mgi/Mg2 


BOS x5 


F 


HB06 


0.8 


H07 


Y 


(Mgi/Mg2)'/2 


BOB 


G 


HB06 


0.8 


H07 


Y 


(Mgi/Mg2)2 


BOB 


H 


HB06 


1.0 


H07 


Y 


Mgi/Mg2 


BOB 


I 


HB06 


2.0 


H07 


Y 


Mgi/Mg2 


BOB 


L 


HB06 


3.5 


H07 


Y 


Mgi/Mg2 


BOB 


Note. — 


Parameter of reference (A) and tost models (B-L). 


Values in bold face 


are differences to the 



reference model. 

Results shown in Figures lA.ll IA.2I and IA.3I 

Redshift dependent SFR; HB06 refers to IHopkins fc BeacomI (200B) 
Exponent of mass dependence in SFR (see Eq. |2) 

Redshift dependent AGN luminosity; H07 refers to IHopkins et al.l I I2007I') . Figure 8 
Stochastic element in BH growth 

Amount of dynamical friction in galaxy mergers; BOS refers to lBovlan-Kolchin et all 1(2009 ) 

Finally, models H-L (Figure IA.3[) test the importance of a dependence of SFR on stellar mass. In these models 
we vary the q exponent in equation [2] models H and I do not show any appreciable variation when compared with 
model A. Model L which has an absolutely unrealistic value for q is the only model where we were able to break the 
Mhuige — Mbh correlation. This result can be easily understand in the following way: if star formation is a too strong 
function of stellar mass, then the vast majority of stars will be formed in the main branch of the tree that hosts (by 
definition) the most massive progenitor of our galaxy. This implies that the bulk of stellar mass will not be subject to 
any averaging process and hence the central limit theorem does not apply. Moreover given the artificially high fraction 
of stars produced within the central galaxy, there will be no major merger, this explains why in model L we do not 
get any bulge more massive than 4 x 10^° Mq. 

All other model are not distinguishable from the reference model A, this underlines one more time the convergence 
power of galaxy mergers and shows that the actual implementation of star formation, dynamical friction, BH accretion 
and bulge formation are only secondary effects. 

SCALING RELATIONS WITH STELLAR MASS AND HALO MASS 

Hierarchical assembly produces correlations between any two parameters that take part in this cascade. The main 
focus of this paper lies on the MBH-Afbuige-relation, but AJb h also correlates with total stellar mass and dark matter 
halo mass. This is shown in the two panels of Figure IA.4| which include the receipes for BHA. The clearest and 
simplest correlation is the one with A/dm (right side), which in the simple framework of this toy model is very close 
to slope=l, with some scatter and no bending. 

Mbh vs. total stellar mass (inch SF receipe), as shown on the left side of Figure [K^ is at the massive end largely 
identical to the MBH-Afbuige-view, because most galaxies there will have had sufficient numbers of minor and major 
mergers in order to make them bulge dominated. At the low-mass end, the scatter is still larger than at high masses, 
due to the smaller number of (averaging) past mergers for each halo, but it is somewhat smaller than for bulge mass, 
since the extra random element from disk-to-bulge conversion is absent. This also explains the missing plume of 
seemingly high MeH-systems at low masses visible in Figure 21 which in this way can be explained to actually be disk 
galaxies with small bulges and normal sized black holes for the amount of total stellar mass. 
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Fig. A.l. — Models A-D; see Table[T]for model parameters. 
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Fig. A.2.— 



Models E— G; see Table [T] for model parameters. The reference model A is shown again as a comparison. 
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Fig. A.3.— 



Models H-L; see Table [T] for model parameters. The reference model A is shown again as a comparison. 
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Fig. a. 4. — Left: Black hole vs. total stellar ma ss for 400 of our model halos (red points), compared as in Figure|3]to the observed local 
relation in black, including the compilation from l IFeoli fc Mancirnl 120091 . circles), low-mass spheroids (open squares) and upper limits for 
spiral bulges (triangles) from Greene ct al. (200^). As before, the solid line is the linear fit by Haring & Rix (2004) with a slope of 1.12. 
Right: Black hole vs. dark matter halo mass. The dashed line with slope=l is plotted to guide the eye. Please note: The abcissa is shifted 
by 3dex with respect to the stellar mass plot on the left. 



